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Fake evolution of dark energy from observation data 
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The equation of state (EOS) of the dark energy is the key parameter to study the nature of the 
dark energy from the observation. Though the dark energy is found to be well consistent with 
the cosmological constant with a constant EOS of —1, weak evidences from different observation 
data and analyses show that dark energy models with an evolving EOS slightly less than —1 at 
some medium redshifts and greater than —1 at high redshifts are more favored. In this paper, It 
is shown that how such a pattern of an evolving dark energy EOS can be just biases arising from 
the statistical method widely adopted in data analyses together with the dependence of the cosmic 
expansion on the dark energy EOS. The issue is actually not limited to dark energy or cosmology. 
It represents a class of mathematical problems of Bayesian analysis. It should be paid attention to 
in similar data analyses to avoid biases in drawing conclusions. 


I. INTRODUCTION 


Since the discovery the accelerating expansion of the 
universe Qi , which is attributed to the dark energy, the 
equation of state (EOS) of the dark energy have been, 
and will be for a long time in the future, the key pa¬ 
rameter to study the nature of the dark energy from the 
observation. Though a constant dark energy EOS of — 1 
which corresponds to the cosmological constant is well 
consistent with the observation (see e.g. ii for recent 
analyses), weak evidences are persistently observed from 
different data sets that the dark energy EOS is slightly 
less than —1 at medium redshifts and greater than — 1 
at high redshifts. These evidences come from indepen¬ 
dent analyses of different kinds of data sets, including 
both standard candles like type la supernovae (SNe la) 
and gamma-ray bursts (GRBs) and standard rulers like 
baryon acoustic oscillations (BAO). See @ for a compre¬ 
hensive study showing this trend with SNe la, BAO, and 
other data sets and [a-Q for a series of studies showing 
that GRBs favor a dark energy EOS greater than — 1 
at high redshifts. See also [1, [l^ for some recent stud¬ 
ies with similar results. Interestingly, even mock data 
generated by assuming a AGDM cosmological model also 
show weak evidences for such a trend [111] . In this pa¬ 
per, It is shown that how the trend can be just biases 
resulted from the statistical method widely adopted in 
data analyses together with the dependence of the cos¬ 
mic expansion on the dark energy EOS. 
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II. LUMINOSITY DISTANCE AND STANDARD 
CANDLES 

We start from constructing the procedure for estimat¬ 
ing constraints of standard candles on cosmological pa¬ 
rameters using mock data. Consider a luminosity relation 
of the following form: 

y = Co + ^ CiX^ + e, (1) 

i 

where XiS are some luminosity indicators which can be di¬ 
rectly measured from observation, e is a random variable 
accounting for the intrinsic scatter aint of the relation, 
and y has the form of 

y = log (iirdlK) , (2) 

where dz, is the luminosity distance and J- may be any 
physical quantity that can be directly measured from 
observation. The luminosity relation of Eq. m has in¬ 
corporated all the GRB correlations summarized in [l^ . 
as well as the relations used to derived distance moduli 
from SNe la. For example, the SALT2 method gives 
piB = rrig — M + a ■ xi — fj ■ c. Let y = (ys — m^)/2.5, 
it can be rewritten in the form of Eq. O- Note that the 
intrinsic scatter in distance modulus should be divided 
by 2.5 to convert to the intrinsic scatter aint in V- 

From a similar derivation like that in [l^ . we know 
that, for a given sample of standard candles, the joint 
likelihood function for the coefficients c, the intrinsic 
scatter (Tint, and the cosmological parameters 9 is 

£(c, (Tint, 9) = k'^Ci, (3) 

i 

where k is the normalization factor, i runs over the stan- 
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dard candles, and 
C,= 




X exp 


{y^ - Co - Ej Cja^j.i) 

(^int + + E, c2a2 . 


( 4 ) 


with j running over the luminosity indicators in Eq. O- 
This likelihood function is a quite general one with the 
well known statistic being a special case of it. 

In the discussions below, we ignore the measurement 
uncertainties ^ and ay^ in Eq. O) to simplify the 
problem. 

First consider the simplest case. For the luminosity 
relation of 


y = a + hx + e, (5) 

to estimate its constraints on the cosmological parame¬ 
ters 0 using mock data, we follow the following steps: 

1. Set the hducial values for the parameters of the 
luminosity relation a, 6, cri„t and the cosmological 
parameters 9. Here, we use ao, bg, (Jint.Oj and Og 
to denote the hducial values for the corresponding 
parameters. 


and use li as the abbreviation for l{zi,9,9g). So, 
we have 


2/i — Oo -f bgXi -\- £i -\- li- 


3. Calculate the likelihood. Ignoring the measurement 
uncertainties, we have 


^ 1 
Li = -exp 

Cint 


Cint 


■ exp 


{yi-a - bxjf 

2ct.^ t 
int 

{5a + 5bxi — £i — lif' 


where 5a = a — ag and 5b = b — bg. Let rrii = 
5bxi — £i — li, the joint likelihood function is 


N 


C{a,b,crint,9) = 


-^exp 

hnt 

k 




exp 


2=1 


T,i=i i5a + m^) 

^^int 


(da)^ -I- 2m5a + 'm? 




where m and im? are the average values of m and 
over the standard candles. 


2. Generate the mock data for a sample of standard 
candles. Assume the total number of standard can¬ 
dles is N. For the ith standard candles, we gener¬ 
ate its redshift Zi and luminosity indicator xt and 
draw a sample £i from the distribution of random 
variable £, i.e., the normal distribution M{Q, g)- 
From them we know its Hducial value for y. 


J/i.o = ao + bgXi + £i. 

Then we calculate iFi, i.e. the value of T for indi¬ 
vidual standard candles that should be measured 
from the observation given the above information, 
through 

2 /i.o = log [Andlizi, . 

Thus, we can calculate yi for any given cosmological 
parameters: 


4. Marginalize out the nuisance parameters. Since we 
aim at constraining the cosmological parameters 9, 
we need integrate the joint likelihood function over 
the nuisance parameters. In this case, they are the 
calibration parameters a, 6, and (Tint- 


(a) Integrate over the intercept parameter. 


/ -I-OO 

C{a,b,aint,9)<15 

-OO 






where cr^ = ^ Efci(^i ~ 1® 

the variance of m. 

(b) Integrate over the slope parameter. Let tii = 
£i + li, then rrii = 5bxi — rii. Since Xi, £i, and 
li are independent of each other, we have 


Vz = log [‘iTrdl{zi,9)J^i] 

= log [4TTdl{Zi,9g)I'i] 

-log [dl{zi,9g)] -flog [dl{zi,9)] 

= yifi -f 21og [dL{zi,9)/dL{zi,9g)] 

_ , ^ , , 01 dL{zi,9) 

— Oq -f bgXi + £i + 2 log ——p - . 

dL(Zi,9g) 

For later convenience, we dehne 

l{z,9,9o) = 21og (6) 


a 


2 

m 


2^2 


{5b) a, 


2 

n ■ 


Therefore 


/ + 00 

C{b,a{nt,9)d5b 

-OO 


’ —OO 

k 27r 1 


N-2 


N a, 


• exp 


2<t/^ 


Note that here is the standard deviation of 
{xi}, not the measurement uncertainty i 
in Eq.®. 
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(c) Integrate over the intrinsic scatter. 


j'+oo 

C{6) = / /:(crint,6')daint 
Jo 


kiT N-3 f N — Z 

-2^~r 


Na^ 

K)- 


(Nal) 


N-3 
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Here cr^ = cr^ + af, and cr^ is an estimation 
of cr?jjQ. So the marginal likelihood of the 
cosmological parameters 9 is given by 

C{9) oc {af^t,o + ^f) ^ ■ 


For the general case of Eq. o, the steps and the 
derivation are basically the same as the above. The 
step 113 can be repeated until all the slope parameters 
are integrated over. It is easy to check that the general¬ 
ized marginal likelihood of the cosmological parameters 
9 is given by 


£(0)cx(aL,o + ^f) (7) 

where p is the number of the calibration parameters 
which include the coefficients c and the intrinsic scatter 
CTint- This also applies to the case of no luminosity indi¬ 
cator, i.e., the case of the luminosity relation y = a + e. 

Thus, at this point, given a luminosity relation, we do 
not need to follow the above steps any more to estimate 
its constraints on cosmological parameters. Instead, we 
can directly calculate the marginal likelihood of the cos¬ 
mological parameters using Eqs. ([6]) and 0 with the 
input of the number of luminosity indicators involved, 
the intrinsic scatter of the luminosity relation, and the 
number of the standard candles and their redshifts. 

Throughout this paper to the end, if not stated explic¬ 
itly otherwise, the flat ACDM with = 0.3 is used as 
the fiducial cosmological model. And since our focus is 
on the dark energy EOS wlz), we fix all other cosmolog¬ 
ical parameters at their fiducial values. Thus l{z,9,9o) 
reduces to l{z,w{z),—1). For the wCDM cosmological 
model where w{z) = w is a constant along the redshift, 
it further reduces to l{z,w, —1). 

Now let us take a close look at l{z,w,—l). From 
Eq. (13), it is easy to check that l{z,w,—l) approaches 
0 when z approaches 0 and, as long as w < 0, l{z, w, —1) 
also approaches 0 when z approaches infinity. So 
Z(z, in, —1) has a maximum or minimum at some red- 
shift where its differential with respect to z is equal to 0. 
See Eig. [T] (left panel) for examples of l{z,w,—l) versus 
z for some values of w. Eor a given w = Wnp, we can 
find the redshift, say Zdataj where dl{z,w,—l)/dz = 0. 
Since af is the variance of l{z,w, —1) along the redshift, 
if we only use standard candles distributing in a narrow 
redshift range around Zdata 7 then ai, as a function of w, 
will show a local minimum at w = Wnp in addition to at 
the fiducial value m = — 1. So, from Eq. ©, we know 


that the marginal likelihood of re, C{w), will have local 
maxima at both w = Wnp and w = —1. See Eig. [5] for an 
illustration of two examples. The fiducial value w = — 1 
is what we want from the constraining, while w = w^p is 
an irrelative non-physical value. A blind analysis with¬ 
out considering the impact of Wnp could lead to biased 
conclusions in the constraining of the dark energy EOS. 
The relation between w^p and Zdata for standard candles 
is plotted in Eig. [T] (right panel). We can see that Wnp in¬ 
creases with the increase of Zdata and crosses the fiducial 
value w = — 1 at Zdata — 1-3. This means that, taking 
into account both the likelihood peaks corresponding to 
the hducial value —1 and the non-physical value w^p, as 
was done implicitly in related analyses, the standard can¬ 
dles with redshifts less/greater than 1.3 would appear to 
favor a dark energy EOS less/greater than —1. This is 
illustrated in Eig. [2] 



FIG. 1. l{z,w, —1) versus z for some values of w and Wnp 
versus Zdata for standard candles. 



FIG. 2. at versus w and the corresponding probability dis¬ 
tribution of w. For the left column, 100 standard candles 
uniformly distributing in the redshift range [0.5,1] are used, 
and (Tint.o = 0.06 and p = 4 are assumed for the luminosity 
relation. For the right column, 100 standard candles uni¬ 
formly distributing in the redshift range [2, 4] are used, and 
Uint.o = 0.1 and p = 3 are assumed for the luminosity relation. 

With that the likelihood peak corresponding to Wnp 
evolves with Zdata while the one corresponding to the 
fiducial value does not, the impact of w^p actually would 
not be an issue if we use a sample of standard candles 
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with a wide redshift distribution and assume a constant 
dark energy EOS, since the likelihood peak corresponding 
to the hducial value would be strengthened by standard 
candles at different redshifts while the one correspond¬ 
ing to Wnp would be suppressed. However, when a flex¬ 
ible parameterization of the dark energy EOS allowing 
evolution along the redshift is used, which is a realistic 
demand to reconstruct the dark energy EOS from obser¬ 
vation, the impact of w^p enters since, for these cases, the 
constraint of the dark energy EOS at a certain redshift 
mainly comes from only part of the standard candles at 
some redshifts. Eor Zdata near to 0, w^p is far below the 
fiducial value w = — 1, its impact in constraining usually 
has been implicitly eliminated by assuming that the dark 
energy EOS is not very far away from —1. Thus, if the 
impact of Wnp is significant enough, we can expect, from 
the constraining, a dark energy EOS that is around — 1 
at low redshifts, slightly biased to below —1 at medium 
redshifts, and slightly biased to greater than —1 at high 
redshifts. This is exactly the trend that is observed in 
different analyses [MU- The significance of the biased 
trend depends on, in addition to the redshift distribu¬ 
tion of the data, the flexibility of the parameterization 
of the dark energy EOS. A less flexible parameterization 
means a stronger prior introduced, which may erase the 
bias if it is strong enough (an extreme example is the pa¬ 
rameterization of a constant dark energy EOS along the 
redshift). However, a too flexible parameterization usu¬ 
ally means large statistical errors, which could totally 
overwhelm the bias and make it not notable. Especially 
for cases with a poorly constrained dark energy EOS at 
high redshifts, the likelihood there usually has a long 
wing at the lower end (since the data does not care much 
about the value of the dark energy EOS at high redshifts 
as long as the matter dominance is guaranteed, which 
means only a upper limit is imposed on the dark energy 
EOS), which could make the dark energy EOS at high 
redshifts appear less than —1 instead of greater than — 1. 
The relation between the significance of the trend and the 
parameterization of the dark energy EOS is very subtle. 
A principal component analysis |15| of 1 -I- w{z) utilizing 
Eqs. dH]) and ([7]) is presented in Eig. [31 It shows that the 
best constrained components do reflect the impact of Wnp 
with w biased toward less than —1 at medium redshifts, 
which is consistent with the result of similar analyses 
derived from real data M- (The best constrained com¬ 
ponents do not include the variation of the dark energy 
EOS at high redshifts, so the bias from Wnp at high red¬ 
shifts does not show up here.) It can be expected that, 
if the prior introduced by some parameterization of the 
dark energy EOS and/or extra data effectively suppress 
the poorly constrained components while retaining the 
best constrained components, the bias from w^p will be 
visible in the final result. Eor example, Eig. 2 (Panels Al, 
A2) in @ shows that only the 2(3) strongest data modes 
survive the addition of the prior there, thus, when com¬ 
pared with the result here derived from the fiducial model 
of the cosmological constant, the favor of the dynamical 


dark energy in Q is understandable as a result of the 
impact of Wnp- 

It should be also noted that, when using Eqs. m and 
the problem has been simplified in a few places. For 
example, the measurement uncertainties are ignored and 
the covariance between x, e, and I are assumed to be ex¬ 
actly zero (Though it is expected to be true, real data are 
usually not so ideal). How these simplifications impact 
on the constraining of the evolution of the dark energy 
EOS need further investigation. For a precise compari¬ 
son with results derived from real data, they may have 
an unignorable effect and should be taken into account. 
And in practice, standard candles are usually handled 
differently in different analyses before used to constrain 
cosmological parameters, which is another thing need to 
be considered. 


III. GENERALIZATION 


It is easy to see that the above discussions can be ap¬ 
plied to more than just luminosity distance and standard 
candles. In fact, Eq. ([3]) can be generalized to 

y = \og[X{z,e)F], ( 8 ) 


where X{z, 9) can be anything to be measured from the 
observation that depends on cosmological parameters. 
All the subsequent derivations are of the same, except 
that the definition of l{z, 9, 9o), Eq. ([S]), is generalized to 


l{z,9,9o) 


log 


X{z,9) 

Xiz,9o)' 


(9) 


The form of Eq. ([7|) remains unchanged. Thus, the results 
are generalized to other observations. 

For BAO survey, the comoving sound horizon at the 
baryon drag epoch, rs(zdrag), is used as the standard 
ruler. For the standard ruler of the transverse direction, 
of the line-of-sight direction, and of the combined direc¬ 
tions, the following distance ratios are measured from the 
survey respectively: 


9s 


Szs 

ds 


^s(.^drag) 

(1 -|- z)dA 

‘I"s i^drag) Hi^z') 
C 

l^si^drag) 

Dv 


( 10 ) 

( 11 ) 

( 12 ) 


These relations can be rewritten in the form oiy = a-\- e 
with a = log[rs(zdrag)], X being 9s, 6zs, and ds respec¬ 
tively, and X{z,9) being (I -|- z)dA, c/H{z), and Dy 
respectively. Thus, l{z,9,9o) can be calculated from 
Eq. (|ni) and corresponding Wnp can be derived from it. 
{Here, only BAO measurements from density fluctua¬ 
tions of baryonic matter are considered, no prior in¬ 
formation about rs(zdrag) from cosmic microwave back¬ 
ground (CMB) measurements is inputted, rs(zdrag) is 
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FIG. 3. Principal component analysis El of 1 + w{z). w{z) was parameterized using 20 bins with a constant bin width of 
0.03 in scale factor within [0.4, 1] and jl + w{z)\ < 50 was assumed. Mock SN la data (crint.o = 0.06 and p = 4) together with 
mock GRB data (crint.o = 0.4 and p — 3) were used. The mock SN la data has the same redshifts of SNLS3 data [Q and the 
mock GRB data has the same redshifts of the sample in [^. From left to right, the 4 best coinstrained components of 1 + w{z) 
are plotted in the first panel (Small shifts in 2 are applied to the components to reduce the overlap of the lines so that the 
components are presented more clearly). In the following panels, the median values (red lines) and the Icr confidence intervals 
(blue lines) of w{z) reconstructed from the 2, 3, and 4 best constrained components of 1 + w(z) are plotted. 


treated as an unknown constant and simply marq inal- 
ized during the constraining. For example, in til . I la/ . 
part of the BAO measurements are summarized in dis¬ 
tance ratio Hy (0.35)/Hy (0.2), which is effectively equiv¬ 
alent to marginalizing out rs{zdra,g) if used to constrain 
the dark energy EOS. If BAO measurements are com¬ 
bined with CMB measurements, such that the informa¬ 
tion about ra{zdra.g) is inputted, the impact ofwnp should 
be relieved or eliminated.) Since the luminosity distance 
dp and the angular diameter distance dA relate to each 
other though dp = ^^(1 + 2 )^, the l{z, 9,9o) for the stan¬ 
dard ruler of the transverse direction differs from that for 
standard candles only by a constant factor of 2. So the 
derived w^p is totally the same. For the other two cases, 
i.e., the standard ruler of the line-of-sight direction and 
of the combined directions, Wnp versus ^data is plotted in 
Fig. m We can see similar behaviors of Wnp versus ^data- 
Wnp crosses the fiducial dark energy EOS at different val¬ 
ues of Zdata- For the former case, at Zdata — 0.65, and for 
the latter one, at Zdata — 0.9. 




FIG. 4. Wnp versus Zdnta. for the standard ruler of the line-of- 
sight direction (left) and of the combined directions (right). 

To be more general, what were discussed here is not 
limited to dark energy or cosmology. In fact, it represents 
a class of mathematical problems of Bayesian analysis. 
Fitting data to a linear relation and using the relation 
to constrain model parameters are very common in data 
analyses. The key quantity in the discussions, ai, appears 
during the marginalization (it shows up in|4^through am 
when we integrate over the intercept parameter), which 


is almost inevitable in data analyses. So similar biases 
may show their shadow here and there in similar prob¬ 
lems. Concerning data analyses in cosmology, a lot of 
parameters are involved. Marginalization is used very 
often. The kind of bias discussed here should be checked 
carefully in the era of precision cosmology. 


IV. SUMMARY 

Starting with luminosity relations of standard candles, 
the steps was described that one should follow to esti¬ 
mate the constraints of a luminosity relation on cosmo¬ 
logical parameters using mock data, from which a simple 
and general formula was deduced that can be used to di¬ 
rectly calculate the marginal likelihood of cosmological 
parameters. Using the formula together with the depen¬ 
dence of the cosmic expansion on the dark energy EOS, 
it was shown how a kind of bias could arise that leads 
to a fake evolution of the dark energy EOS, whose sig¬ 
nificance depends on the flexibility of the parameteriza¬ 
tion of the dark energy EOS and the redshift distribution 
of the data. Then the formula was generalized to more 
than just standard candles. It was shown that the BAO 
data could lead to similar biases. It was mentioned that 
the issue represents a class of mathematical problems of 
Bayesian analysis and should be paid attention to in sim¬ 
ilar analyses. 
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